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Abstract 

The linear programming method is applied to the space U„(C) of uni- 
tary matrices in order to obtain bounds for codes relative to the diversity 
sum and the diversity product. Theoretical and numerical results improv- 
ing previously known bounds are derived. 



1 Introduction 

Nowadays breakthrough of wireless communications has provided new and nice 
problems to the field of coding theory. Indeed strategic issues of MIMO commu- 
nications has lead to consider coding in Grassmann and Stiefel manifolds, and 
also in unitary groups [13j . Codes in unitary groups are useful in the context 
of non-coherent flat Rayleigh channel as shown in [BJ. The performance of a 
unitary space-time code V is measured (see [7]) by two functions namely the 
diversity sum SV and the diversity product IIV : 

£V:= J_niin{||a:-y|| : x ^ y e V} , 
2V« 



-min||det(a;-t/)|« : x ^ y € vj 



IIV := — min 

2 

Here ||A|| denotes the standard Euclidean norm of complex matrices: ||A|| 2 = 
Trace(AA*) = \Ai,j\ 2 - A standard problem is, given a number N of points 
in U„(C), to maximize the value of SV or of IIV. Many authors have addressed 
this question, narrowing gaps between bounds and explicit constructions (see 
0, 0) [ID]) [II])- The linear programming method, which was initially devel- 
oped by Philippe Delsarte in the framework of association schemes [I] , and is a 
powerful method to deal with such questions, has not been applied to unitary 
codes before. Delsarte method was successfully adapted to the compact two- 
point homogeneous spaces by Kabatiansky and Levenshtein ,9| and recently to 
more general situations like the Grassmann codes [TJ, the permutation codes 
[T^] , the ordered codes [5] , [5] . Most of the situations mentioned above fit into a 
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common framework, namely a compact group G acts homogeneously on the un- 
derlying space X, and the representation theory of G constructs a certain family 
of orthogonal polynomials naturally attached to X . Standard methods of har- 
monic analysis show that these polynomials hold the desired positivity property 
that allows for Delsarte linear programming method. This general framework 
is recalled in Section II, and we show in Section III that unitary codes can be 
treated likewise, the Schur polynomials being the associated family of orthog- 
onal polynomials. Section IV and V present the results, both numerical and 
analytic, obtained by the implementation of this method. It turns out that we 
improve all previously known bounds concerning the diversity sum and the di- 
versity product. Moreover it is worth pointing out that the mathod can easily 
be extended to more complex situations, for example a diversity function in- 
volving both YjV and IIV. 



2 The linear programming method on homoge- 
neous spaces 

We briefly describe the linear programming method on homogeneous spaces. 
For more details we refer to [9], [HI Chapter 9] for a treatment of 2-point 
homogeneous spaces, and to [T] for the prominent case of Grassmann codes. 

Let G be a compact group acting transitively and continuously on a compact 
space X, and r : X x X — > Y such that r characterizes the orbits of G acting 
on X x X. We mean here that, for all x\, X2,x[, x' 2 E X, 

t(x 1 ,x 2 ) = t(x{,x' 2 ) 3g e G : g(xi,x 2 ) = {x[,x' 2 ). 

Let S be any subset of Y, we call a finite subset V C X a S-code if for all 
ci ^ c 2 € V, t(ci,c 2 ) € S. 

A continuous function P : Y — > C is said to possess the positivity property 
if for any finite subset V C X and any complex function a : X — > C, 

a(x)a(y)P (t(x, y)) is real non-negative. 

x,y£V 

A canonical example is the constant function Pq = 1; non-trivial examples are 
given by the so called zonal functions that we introduce now. 

Let L 2 (X) — {u : X ~ * C : J x \u(x)\ 2 dx < oo} where dx is the unique G- 
invariant Haar measure on X such that J x dx = 1 . This vector space is given 
the standard G-action defined by g.uix) = u(g^ 1 (x)) and is endowed with the 
canonical G-invariant hermitian product : (ui,u 2 ) = j x u 1 {x)u 2 {x)dx . The 
Peter- Weil theorem shows that L 2 (X) can be decomposed as a direct sum of 
G-irreducible subspaces Vj. The next step is to associate to each irreducible 
subspace Vi a so-called zonal function Py i . A standard construction is the 
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following: given an orthonormal basis (ui(x), ...,Ud(x)) of Vi, one can define 



1 p 

Pvi(x,y) = - 22ui(x)ui(y). (1) 

di — f 

2 — 1 

where di = dim(Vi). Since these functions are constant on G-orbits we can 
rewrite Pvi{x,y) = Py i (r(x,y)). From this property comes the term zonal 
functions used to qualify them. From equation ([1} it is easy to prove that these 
zonal functions verify the positivity property and do not depend of the chosen 
orthonormal basis. It turns out that, when the irreducible subspaces Vi are 
pairwise non isomorphic, the cone of continuous positive G-invariant functions 
is exactly the set of linear combinations with non negative coefficients of the 
Py i (see [S]). In the remaining of this paper we assume that this condition is 
satisfied. We moreover let Vb denote the one-dimensional subspace associated 
to the trivial representation of G. 

The so-called linear programming bounds are obtained with the following 
theorem : 

Theorem 2.1 Let P = CjPvi a linear combination of the zonal functions 
Py i with a finite number of non zero coefficients. Assume furthermore that : 
Cj > 0, Co > and 5R(P) (the real part of P) is non-positive on S. Then any 
S-code verifies 

\V\ < ^) ( 2 ) 
Co 

where tq ~ r(x, x) for any x. 
Proof : On one hand, 

J2 P(r(x,y))= J2 nr{x,y))+ £ P(r(x,y)) 

x,y£V x—y^V x^y£V 

< |V|P(r ) 

On the other hand, 

P(r(x,y)) = coPo(r(a;, y)) + non negative terms 

x ,y£V x ,y£V 

> |v|V 



3 The case of unitary codes 

As a particular case we set G = U„ (C) x U„ (C) and X = U„ (C) . For g = U x V 
and x = M we set gx — UMV^ 1 . In this context the orbit of a pair (x, y) is 
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characterized by the eigenvalues (e* 6 * 1 , e lSn ) of the unitary matrix xy 1 . The 
Peter- Weil theorem gives a decomposition of L 2 (X) into irreducible subspaces 

L 2 {X)=@{V X ®V X ) (3) 

where the sum runs over all irreducible representations V x of U„ (C) . It is worth 
noticing that the G-subspaces V x ® V x are pairwise non isomorphic. From this 
decomposition one can deduce the following theorem : 

Proposition 3.1 The zonal functions associated to this decomposition are 

where x denotes any irreducible character of U„(C). 

The irreducible characters of U„(C) are known to be finite dimensional and 
to have a nice description using Schur polynomials (|16j). We recall briefly some 
notations and definitions concerning those polynomials. 

For any integer k, a partition of k in n parts is a finite decreasing sequence of 
n non negative integers which sum exactly to k (k is also called the degree of the 
partition). A partition A is dominating (noted A >- n) if Vr < n, Y^i=i ^» — 
Si=i Mi- Given a partition A = [Ai, A„] we define the elementary symmetric 
polynomials m\ 6 Z[xi, ...,x n ] as the renormalization to a monic polynomial of 

Ei, A 2 A„ 
X cr(l) X a(2)--- X a(n)- 

ff€6„ 

These polynomials form a basis of the set of symmetric polynomials Z[xi, x n ] Sym . 
Schur polynomials have been intensively studied and have several definitions. 
For our purpose we will define the Schur polynomials S\ as 

S\ = ^ 

where the Kx ^ 6 N are the so called Kostka numbers. For more precision on 
those numbers see |15j . It is clear that Schur polynomials form another basis of 

z[ Xl ,...,x n ] s y m . 

It is well known that the irreducible polynomial characters of U„(C) are 
expressed using Schur polynomials (we refer to [16j for details) in the following 
way: let (e 101 , ...,e l8n ) denote the eigenvalues of a unitary matrix M, and A a 
partition. Then 

Xx (M) = S x (e ie \...,e w "). 

One obtains all irreducible characters of U„ (C) by multiplying the characters x\ 
by a relative power of det(Af) = Y[ e l9k . All together, we obtain the theorem: 

Theorem 3.2 For all partition A and s£Z, let 

P\, S {XI, ...,X n ) = (xi., .XnYPxixx, . . .,X n ). (4) 
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These rational fractions give the zonal functions associated to the irreducible 
decomposition ([3]) in the following way: if x — det s <8>A, if the eigenvalues of 
xy^ 1 are (e 101 , e l8n ), then 



We are now almost ready to compute bounds for unitary codes V C U„(C). 

Let 

1 ™ 

dUe l9 \...,e^):=-J2(l-co S 9<) 

i—l 

and 

1 / - 

\i=l 

and define a%(x,y) (resp. d^ix, y)) to be the above functions evaluated at 
the eigenvalues of xy^ 1 . These functions are related to the diversity sum and 
diversity product by 

SV = min ds(x,y) and LTV = min du(x,y). 

x^y x^y 

We recall that the orbit of a pair (x, y) is characterized by the eigenvalues of 
xy -1 , so ds and du are constant on G-orbits. We may now define, with the 
notations of Section II, the sets S related to each diversity function: 

S s (8) :={(e i0 \...,e ie »):d s (e i0 \...,e i9 «)>5} 

Sn(S) := {(e i( \..., e i9 ») : d u (e i6 \ e ie ") > 6}. 




4 Analytic bounds of low degree 

From the explicit description of the zonal functions (H]), we have deduced con- 
venient polynomials of low degree which verify the positivity property. Using 
formula we derive the analytic bounds of Theorems 14.11 and 14.21 



Theorem 4.1 Let V be a unitary space time code with diversity sum SV, the 
following upper bounds hold : 
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IVI < 



16n 2 (£V) 2 



2n(2n 
*/(£V) 2 



1)(EV) 2 - 
2n 2 



{2n 2 



(7) 



> 



2n(2n - 1) 



Theorem 4.2 Let V be a unitary space time code with diversity product T1V , 
the following upper bounds hold : 



8n(nV) 2 , n9 2n-l 

' v ' * M nv) 2 -(2n-i) ' 4/ (HV) > — ' ^ 3 (9) 

8(nv) 6 + 4(nv) 4 + 8(nv) 2 
8(nv) 



| V | < V ^ o^6_T . (io) 



«/(nV) 2 > i n = 2. 



The proofs of these theorems are based on Theorem l2.1l and on the following 
lemma : 

Lemma 4.3 Let yj = cos 6^- and rnu lt . ^ ar \(y) the elementary symmetric poly- 
nomials in the yj — cos 9j . The following polynomials are linear combination of 
the zonal functions ((4]) with non negative coefficients: 

Q[o] = i 

Q[u] = m [n]{y) + 

Q[i] = m {1] {y) 

Q[2) = m [2] (y) +m [u] {y) - 

Q[2,i] = m [2,i] (y) + 2m [i4,i] (y) - i m [i] (y) 

Q[3] = m [3]{y) +TO[ 2 ,i](y) +m [liljl ](y) - ^^pm^y) 



Proof : [Theorem 01 Sketch] Let s = ~ £™ = i cosfl,- so that s = 1 - 2d|. We 
apply formula |2]) with the polynomials : P = Q^—ns, P = (Qm— ns)(Qm+n), 
P = (Qm — ns)R, where i? is the symmetrization of (yi + l)(j/2 + 1). 



G 



Sketch] Let p = 2d*. 
Q m - n(p - 1), P = Q [2] 



Proof : [Theorem 
polynomials : P = 
(if n > 3), 

P = Q[2] + {^ + 2p- l)Q m + (p 3 - |) (if n = 2). 



We apply formula ^ with the 

(n+l)p n (n+l)(2n(p-l) + l) 
H o Vl H 4 



5 Numerical bounds 

Numerical programs give accurate approximations of the best linear program- 
ming bounds over a large interval of validity, not covered by the bounds proved 
in Section IV. The following curves plot the linear programming bound on the 
diversity functions as a function of the cardinality of the code. The programs 
optimize the choice of a polynomial in the variables (cos#i, . . . ,cos#„), with 
degree at most equal to some parameter D. Increasing D gives accurate re- 
sults over a wider range of values for the diversity functions, but also increase 
the computational time. We use D = 19 for n — 2 and D = 13 for n = 3. 



Bounds for the Diversity Sum, n-2 
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Bounds for the Diversity Product, n=2 



Numerical degree=7 
Numerical degree=19 ■ 
Bound from d y (deg=1 9) - 



100 
Cardinality 

Bounds for the Diversity Sum, n=3 



Num. degree=7 
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N 


24 


48 


64 


80 


bounds in [TU] 


0.6746 


0.6193 


0.5969 


0.5799 


Bl 


0.7598 


0.6603 


0.6131 


0.5932 


B2 [7] 


0.7794 


0.6734 


0.6235 


0.6026 


LP cfe 


0.6547 


0.5797 


0.5488 


0.5254 


LP d n 


0.5730 


0.4989 


0.4711 


0.4504 




N 


100 


120 


128 


1000 


bounds in 10J 


0.5632 


0.5499 


0.5452 




Bl [7] 


0.5578 


0.5425 


0.5347 


0.3270 


B2 [7] 


0.5654 


0.5496 


0.5415 


0.3285 


LP cfe 


0.4999 


0.4816 


0.4753 


0.2964 


LP d n 


0.4301 


0.4144 


0.4089 


0.2574 



Table 1: n = 2 



N 


24 


48 


64 


80 


LP d s 


0.7178 


0.6939 


0.6797 


0.6692 


LP d n 


0.6431 


0.5942 


0.5752 


0.5628 




N 


100 


120 


128 


1000 


LP <fa 


0.6598 


0.6532 


0.6511 


0.5586 


LP d n 


0.5482 


0.5369 


0.5332 


0.4330 



Table 2: n = 3 



The Tables I and II compare the linear programming bounds for the dimen- 
sions 2 and 3 with the previous results of |10j and [7] and show an improvement 
in all cases. 

These tables give upper bounds for the diversity when the cardinality N is 
fixed. All entries except the ones of the last line (LP dn) are concerned with the 
diversity sum. The second line tabulates the bounds settled in [10], obtained 
using Coxeter upper bounds. The two bounds of [7] were obtained using sphere 
volume computations. 

Moreover, concerning the diversity product, both our numerical results and 
analytic results (compare ([6]) and ([9])) show a large gap between the bounds 
for diversity sum and diversity product, in favor of the diversity product. This 
is worth to point out since in previous publications bounds for the diversity 
product were essentially deduced from the trivial inequality SV > IIV and 
hence appear to be weak. 

6 Conclusions 

In this paper we have developed the linear programming method for the unitary 
space time codes. We have obtained both numerical and analytic bounds. The 



9 



results improve previously known bounds. Furthermore the linear programming 
method allows to deal with non-distance functions as the diversity product di- 
rectly. 
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